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ABSTRACT 

Swift discovered XRF 050416A with the Burst Alert Telescope and began observing it with its 
narrow field instruments only 64.5 s after the burst onset. Its very soft spectrum classifies this event 
as an X-Ray Flash. The afterglow X-ray emission was monitored up to 74 days after the burst. The 
X-ray light curve initially decays very fast (decay slope a ~ 2.4), subsequently flattens (a ^ 0.44) and 
eventually steepens again [a ^ 0.88), similar to many X-ray afterglows. The first and second phases 
end '^172 and ~ 1450 s after the burst onset, respectively. We find evidence of spectral evolution from 
a softer emission with photon index F ~ 3.0 during the initial steep decay, to a harder emission with 
r ~ 2.0 during the following evolutionary phases. The spectra show intrinsic absorption in the host 
galaxy with column density of ~ 6.8 x 10^^ cm~^. The consistency of the initial photon index with the 
high energy BAT photon index suggests that the initial fast decaying phase of the X-ray light curve 
may be the low-energy tail of the prompt emission. The lack of jet break signatures in the X-ray 
afterglow light curve is not consistent with empirical relations between the source rest-frame peak 
energy and the coUimation-corrected energy of the burst. The standard uniform jet model can give a 
possible description of the XRF 050416A X-ray afterglow for an opening angle larger than a few tens of 
degrees, although numerical simulations show that the late time decay is slightly flatter than expected 
from on-axis viewing of a uniform jet. A structured Gaussian-type jet model with uniform Lorentz 
factor distribution and viewing angle outside the Gaussian core is another possibility, although a full 
agreement with data is not achieved with the numerical models explored. 
Subject headings: gamma rays: bursts; X-rays: individual (XRF 050416A) 



1. INTRODUCTION 

, Th e Swift Gamma-Ray Burst Explorer l)Gehrels et al.l 
' 12004(1 ■ successfully launched on 2004 November 20, is 
\ dedicated to the discovery and study of gamma-ray 
bursts (GRBs) and their X-ray and optical afterglows. 
Its fast-pointing capability, compared with previous 
satellites, allows us to repoint towards GRB sources ap- 
proximately 100 s after the burst detection bv t he Burst 
' Alert Telescope fBAT: IBarthelmv et al.ll2005a|) and to 
study, for the first time, the early phases of the after- 
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glow evolution. Moreover, the very broad energy cover- 
age allows a simultaneous study of the phenomenon in 
the optical, soft and hard X-ray bands. 

One of the main Swift results has been the direct 
observation of the transition between the prompt and 
the afterglow emission. A growing number of early 
rapidly-fading X-ray afterglows have been observed by 
Swift, that have been successfully in terpreted as the 
tail of the prompt GRB emission (e.g. Taeliafcrri et aU 
p)05; Cusumano et al. 2006 a; Bart hclmv et al. 2 005M 
Nous ek"et al.l T2006: Chincar ini et all , 2005: Zhang et alJ 
l2006iL The emergence of the true X-ray afterglow com- 
ponent has been identified with a break in the X-ray light 
curve to a less steep (often flat) decay rate, sometimes 
accompanied by spectral variation across the break. Fur- 
ther breaks in the X-ray light curve have been related 
to standard afterglow evolution. The light curve of the 
X-ray counterpart of XRF 050416A show all these char- 
acteristic features. 

BAT detected and located XRF 050416A on 2005 
April 16, 11:04:44.5 UT, at the coordinates RAj2ooo = 
12i^33™5?;6, Dec.T^nnn = -H21°03'10'.' 8, with an uncer- 
tainty of 3' l|Sakamoto et alJl2005alOl . The light curve 
showed a single peak followed by a small bump with 
duration r9o=2.4±0.2 s, with most of the energy emit- 
ted in the 15 — 50 keV band. The time-averaged en- 
ergy distribution was well described by a power law 
{N{E) oc E-'^) with photon index F = 3.1 ± 0.2 (90% 
confidence level) l|Sakamoto et alJl2005(j) . The soft spec- 
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trum and the fact that the fluence in the 15 — 30 keV en- 
ergy band is 6.1 x 10"^ erg cm"^, larger than the fluence 
in the 30 - 400 keV band (1.3 x 10" ^ erg cm"^) clas- 
sify this event as a n X-ray flash (XRF: lHeise 'eraLll2?mt 
iLamb et al. 2005>). A complet e ground analysis of the 
BAT data is presented in .Sakamoto et alJ ((2005c) . These 
authors found that the best flt for the average energy 
distribution of the burst over the T gn interval is given 
by a Band model l|Band et alJll993j) . with peak energy 
Ep = 15.6^2 7 keV, low-energy spectral slope asand ~ 
— 1 (fixed) and high-energy slope /^Band < —3.4 (68% con- 
fidence level). This represents a 3.1 ct improvement with 
respect to a simple power law fit (F = 3.1 ib 0.2 over the 
14—150 keV energy band). Sakamot o et al.l (j2005c) also 
showed that spectral hard-to-soft evolution was present 
during the BAT observation, with the spectrum becom- 
ing considerably softer at the end of each peak, and es- 
timated an isotropic energy Eiso 1.2 x 10^^ erg. 

Following the burst detection, the satellite executed 
an immediate slew and promptly began collecting data 
at 11:05:49 UT (64.5 s after the trigger) with th e Ultravi- 
olet/Optical Telescope CUVOT: lRoming et al.ll2005ft and 
at 11:06:00.6 UT (i.e. 76.1 s after the trigger) with the 
X-Ray Telescope (XRT; Burrows et al. 2005a). 

In the first 100 s of observation UVOT revealed a 
new source in the V filter at RAjaooo = 12'i33'°54!56, 
Dec.T2 nnn = -l-21°03^ 27"3 (with an uncertainty radius of 
0'.'56 (HoUand et al."2006')), with magnitude V = 19.38 
((Schadv ct al. 2005a). Starting with data taken 207 s 
after the trigger, the source was also detected in the U 
and B bands, with magnitudes U — 19.34 ± 0.20 and 
B = 19.85 ± 0.20 mag. There was no further detec- 
tion in the V band down to the 5 a limiting magnitude 
of 19.57 mag (Schadv ct al. 2005b), 173 s after the first 
y-band image. A fading source was also detected at 
193 nm (UVW 2 filter), p lacing an upper limit of 1 to the 
GRB redshift l)Foxll2005j) . The results an d implications 
of the UVOT observations are discussed in lHolland et al.l 
|2006). 

Ground-based optical, NIR, and radio follow-up ob- 
servations were performed with several instruments. A 
fading source w as detected with the A NU 2.3 m telescope 
in the R band (|Anderson et al.ll2005j) . with the Palomar 
200-inch Hale Teles cope in the K^^ ban d ~ 6 min af- 
ter the BAT trigger l|Cenko et al.ll2005ar . with the MAO 
telescope in the R band (i? = 20. 85 ±0.12. with a 900 s 
exposure, 11 hr after the trigger: Ikahharov et al.ll2005() . 
and, marginally, with the KAIT in a 60 s /-band image, 
7.4 min after the trigger (Li et al. 200^. A late observa- 
tion performed with the MAGNUM telescope equipped 
with the MIP dual-beam optical-NIR imager detected 
the afterglow with R = 21.3 mag, 12.2 hr after the 
trigger. Spectra of the host galaxy of XRF 050416A 
were taken with the Low Resolution Imagin g Spectrom- 
eter mounted on the 10-m Keck I telescope ( Genko et al.l 
l2005b'l . The spectrum indicates that the host galaxy 
is faint and blue with a large amount of ongoing star 
formation. Spectral analysis revealed several emission 
lines including [Oil], II/3, H7, and H^, at a redshift 
z = 0.65 35 ± 0.000 2. This is consistent with the pre- 
diction of lFoxl ()2005f l based on the afterglow detection in 
the Swift UVOT UVW2 filter. XRF 050416A is thus one 
of the closest long GRBs discovered by Swift. 



At radio frequencies, no source was detected with the 
VLA down to a limiti ng fiux of 260 uJy at 8 .46 GHz 
37 min after the burst ijFrail fc Soderberdl2005 ^ or with 
the Giant Meter-wave Radio Telescope at 1280 MHz 
~ 9 d after th e burst (p lacing an upper limit of 94 /iJy; 
Ushwara-Ghandra et aljf2005,) . while a source with fiux 
density 260 ± 55 /LtJy was det ected with the V LA at 
4.86 GHz, 5.6 d after the burst llSoderbergll200l . 

In the following, we report on the analysis of the 
prompt emission and of the X-ray afterglow observed by 
Swift. Details on the follow-up XRT observations and 
the XRT data reduction are described in § |21 the tem- 
poral and spectral analysis results are reported in § |21 
In §0]we discuss our results. Gonclusions are drawn in 
§ |S| Finally, in the Appendix we will show how our in- 
terpretation of the early XRT light curve as the tail of 
the prompt emission could be reconciled with the report 
of a peak energy Ep = 15.6^2 7 keV bv iSakamoto et all 
(2005c) and the observed hardness evolution of the BAT 
light curve. 

Throughout this paper the quoted uncertainties are 
given at the 90% confidence level for one interesting pa- 
rameter, unless otherwise specified. We also adopt the 
notation F{y^t) oc t~°'u~^ for the afterglow monochro- 
matic fiux as a function of time, with u representing the 
frequency of the observed radiation and with the en- 
ergy index /3 related to the photon index F according 
to/3 = F-l. 

2. XRT OBSERVATIONS AND DATA REDUCTION 

The Swift XRT is designed to perform automated ob- 
servations of newly discovered bursts in the 0.2 — 10 keV 
band. Four different read-out modes have been imple- 
mented, each dependent on the count rate of the observed 
sky region. The transition betwe en two m odes is auto- 
matically performed on board fsee lHill et a l. 2004. '2005al 
for a detailed description of XRT observing modes). 

XRT was on target 76.1 s after the BAT trigger. It was 
operating in auto state and went through the standard 
sequence of observing modes, slewing to the GRB field 
of view in Low Rate Photodiode (LR) mode, taking a 
2.5 s frame in Image (IM) mode followed by a LR frame 
(1.3 s), and 8 Windowed Timing (WT) mode frames 
(9.6 s), and then correctly switching to Photon Gounting 
(PG) mode for the rest of the orbit. XRT was not able 
to automatically detect the source centroid on-board be- 
cause of its low intensity, but ground analysis revealed 
a fading object identified as the X-ra y counterpart of 
XRF 050416A (IGusumano et al.ll2005^ . XRF 050416A 
was then observed intermittently over 29 consecutive or- 
bits for a total exposure time of 57454 s. During the 8th 
orbit (starting ^ 35 ks after the trigger), the brighten- 
ing of a column of flickering pixels caused uncontrolled 
mode-switching between WT and PG modes. The WT 
data from this orbit are not usable because of their very 
low signal to noise (S/N) ratio. XRF 050416 A was fur- 
ther observed several times up to 74 d later in PG mode. 
The observation log is presented in Table 

XRT data were downloaded from the Swift Data Gen- 
ter at NASA/GSFG (level 1 data products). They were 
then calibrated, filtered, and screened using the XRTDAS 
software package (v. 2. 3) developed at the ASI Science 
Data Genter (ASDG) to produce cleaned photon list 
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TABLE 1 

XRT OBSERVATION LOG OF XRF 050416A 





Sequence 










Exposure 








(yyyy 




1 - d a 11 11 . iTi 111 : s s . s j 


(s since trigger) 




1 


00114753000 


LR'' 


2005- 


-04- 


■16 11:05:49.0 


64.5 


8.2 


1 


00114753000 


IM 


2005- 


-04- 


■16 11:06:00.6 


76.1 


2.5 


1 


00114753000 


WT 


2005- 


-04- 


■16 11:06:08.6 


84.1 


9.6 


1 


00114753000 


PC 


2005- 


-04- 


■16 11:06:18.2 


93.8 


57360 


2 


00114753001 


PC 


2005- 


-04- 


■18 14:27:56.9 


184992.4 


35084 


3 


00114753003 


PC 


2005- 


-04- 


■26 00:53:32.7 


827328.2 


21279 


4 


00114753004 


PC 


2005- 


-04- 


■28 01:07:03.8 


1000939.3 


20710 


5 


00114753005 


PC 


2005- 


■05- 


■02 00:24:42.2 


1343997.7 


6892 


6 


00114753006 


PC 


2005- 


■05- 


■03 00:33:40.8 


1430936.3 


5696 


7 


00114753008 


PC 


2005- 


■05- 


■08 16:38:29.5 


1920825.0 


16702 


8 


00114753009 


PC 


2005- 


■05- 


■13 01:08:16.3 


2297011.8 


22628 


9 


00114753010 


PC 


2005- 


■05- 


■14 01:13:22.3 


2383717.8 


29625 


10 


00114753011 


PC 


2005- 


■05- 


■25 04:03:25.2 


3344320.7 


23897 


11 


00114753012 


PC 


2005- 


■05- 


■26 02:50:04.0 


3426320.5 


692 


12 


00114753013 


PC 


2005- 


■05- 


■27 01:06:28.8 


3506504.3 


5077 


13 


00114753014 


PC 


2005- 


■05- 


■29 01:18:21.3 


3680016.8 


19254 


14 


00114753018 


PC 


2005- 


■06- 


■21 00:46:21.1 


5665296.6 


27979 


15 


00114753019 


PC 


2005- 


■06- 


■22 00:46:54.2 


5751729.7 


24071 


16 


00114753020 


PC 


2005- 


■06- 


■23 00:55:53.8 


5838669.3 


24073 


17 


00114753021 


PC 


2005- 


■06- 


■25 01:07:54.8 


6012190.3 


14581 


18 


00114753022 


PC 


2005- 


■06- 


■28 00:04:59.1 


6267614.6 


13215 


19 


00114753023 


PC 


2005- 


■06- 


■29 00:05:00.2 


6354015.7 


11435 



^The LR observation refers to the settling data set acquired during the slew to the BAT 
coordinates. Only the time interval in which the X-ray countrpart of XRF 050416A is 
clearly detected is considered. 



files^^. The temperature of the CCD was acceptably 
below —50 °C for the whole observation set. The to- 
tal exposure times after all the cleaning procedures were 
8.2 s, 9.6 s, and 368 815 s for data accumulated in LR, 
WT, and PC mode, respectively. 

For both the spectral and timing analyses we used stan- 
dard grade selections: — 12 for PC mode, — 5 for LR 
mode, and — 2 for WT mode. However, data in WT 
mode had insufficient statistics to allow for detailed spec- 
tral modeling and were used only in the light curve anal- 
ysis. Ancillary response files for PC and LR spectra were 
generated through the standard xrtmkarf task (vO.5.1) 
using the response files swxpc0tol2_20010101v007 . rmf 
and swxpd0to5_20010101v007.rmf from CALDB (2006- 
01-04 release). 

In the timing analysis, XRT times are referred to the 
XRF 050416A BAT trigger time T = 2005 Apr 16.461626 
UT (2005 Apr 16, 11:04:44.5 UT). 

3. XRT DATA ANALYSIS 

3.1. SPATIAL ANALYSIS 

Figure d (left panel) shows the XRT image accumu- 
lated in PC mode with a 0.2 — 10 keV energy selection 
during the first and second observations, together with 
the BAT and XRT error circles. The central portion 
of this field is expanded in the right panel of Fig. Q in 
which we show the cumulative image of the follow-up ob- 
servations # 3 to 13. Two sources are visible within the 
BAT error circle. The brighter one is coincident with the 
position of the optical counterpart as derived by UVOT 
(cross point) and it is clearly fading with time. There- 
fore we identify it as the afterglow of XRF 050416A. The 

http: //swift .gsfc .nasa. gov/docs/swift/analysis/xrt_ 
S¥guid6_vl_2 . pdf 



faint er so urce lies 21'.'6 away from the UVOT afterglow 
(|Schadv et a l. 20053,) and it does not show any signifi- 
cant evidence of intensity variations during the XRT ob- 
servations. Its count rate, as determined from the sum 
of the follow-up observations # 3 to 13 (chosen to min- 
imize the contamination from the afterglow emission) is 
(3.3 ±0.7) X 10^4 count s^^. 

The afterglow position derived with xrtcentroid 
(vO.2.7) is RAj2ooo = 12h33"^54!62, Decj2ooo - 
-h21°03'27'.'7, with an uncertainty of 3'.'3. This po- 
sition takes into account the correction for the mis- 
alignment between the telescop e and the satellite op- 
tical axis ijMoretti et alJ I2006|i . The XRT boresight 
corrected coordinates a re 45'.'0 from the BAT position 
iSakamoto et alJl2005bfl and O'.'Q from the optical coun- 
terpart l|Schadv et alJl2005a|) . 

3.2. TIMING ANALYSIS 

The X-ray emission from the counterpart of 
XRF 050416A was detected for the first time in 
the settling data (i.e. those data collected during the 
satellite slew, when the XRT pointing direction was 
less than 10' off the target position). We have a total 
of 8.2 s of data in which the source is significantly 
detected before the "official" beginning of the pointed 
XRT observation. In these data the maximum offset of 
the source is lower than 3' and no vignetting correction 
was thus required. As LR is a non-imaging mode, the 
background has been extracted from pointed LR data 
taken at the beginning of an orbit, 1.5 x 10^ s after the 
burst trigger, when, as we know from the PC data (see 
below), the afterglow emission had faded to less than 
1% of the initial value. This background subtraction 
allows us to correctly account for the emission of the 
serendipitous non-transient sources in the field of view 
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Fig. 1. — Left: Photon Counting mode XRT image of the first and second observations (92 ks exposure), smoothed with a Gaussian 
filter with a 3 pixel standard deviation, showing the XRT and BAT error circles (Moretti ct al. 2006, Sakamoto ct al. 2005Q, respectively). 
The dashed box identifies the area enlarged in the right panel. The UVOT afterglow position is not marked on this image because it is 
indistinguishable from the XRT one. Right: Cumulative im age of the follow-up observations from 3 to 13 (152 ks exposure; Table [Tl . 
showing the XRT error circle, the UVOT position as given in lSchadv et alJ I2005al) (cross) and the position of the second source detected 
within the BAT error circle (square). 



in addition to the instrumental and cosmic X-ray 
background. 

After settling, a single exposure (2.5 s) image mode 
frame was taken, which officially marks the beginning of 
pointed observation. We determined the total amount 
of charge above the background inside a circular region 
of 30 pixel radius centered at the source position in raw 
telemetry units (i.e. Data Number or DN units). We 
obtained 278 DN above a background of 0.5 DN. The 
DN value was converted to count rate by evaluating the 
mean energy of the LR spectrum in the 0.2 — 10 keV 
band (1490 eV). Given a mean energy per DN of ~79 eV 
in the low gain imaging mode we calculated a rate of 
5.9 ± 1.5 count s~^. 

WT data were extracted in a rectangular region 40 
pixels wide along the image strip, which includes about 
98% of the point spread function (PSF). The background 
level was extracted from a rectangular region of the same 
extension, far from the source and affected by minimal 
contamination from other sources in the field. The rapid 
decay of the source and the high background level in 
WT mode allow us to have a significant detection of the 
XRF 050416A only in the first 8 frames (9.6 s) of the 
WT data set. 

Given the pile-up effect in the first part of the observa- 
tion, the presence of the second source and the weakness 
of the afterglow after the first follow-up observation, the 
PC data require different extraction regions for different 
rate levels, in order to optimize the S/N ratio. 

First orbit. The intensity of the source during the first 
orbit of the first observation (Table ^ was high enough 
to cause pile-up in the PC frames. In order to correct for 
this effect, we extracted counts from an annular region 
with an outer radius of 30 pixels (70'.'8) and an inner 
radius of 3 pixels (7'.'1). Such a region includes about 
52% of the PSF. The optimal inner radius was evalu- 
ated by comparing the analytical PSF with the profile 
extracted in the first 2000 s of observation. The ex- 



cluded region corresponds to pixels deviating more than 
one sigma from the best fit of the differential PSF wings 
(i.e. the fit performed on data from 7 pixels outwards). 
The light curve was corrected for the PSF fraction loss. 

First observation (excluding the first orbit) and sec- 
ond observation. In the following 28 orbits of the first 
observation and throughout the second observation the 
intensity of the afterglow was lower than 0.1 count s~^ 
and the pile-up was negligible. The data were extracted 
from the entire circular region of 30 pixels radius in or- 
der to have the maximum available statistics, particu- 
larly important in the last part of the afterglow decay. 
This circle encloses 93% of the PSF. Both here and in 
the previous case, the contribution of the second, fainter 
source detected within the BAT error box is negligible 
with respect to the afterglow intensity. 

Observations 3 to 13. The afterglow had faded to 
a count rate comparable to that of the serendipitous 
source. In order to avoid significant (i.e. greater than 
10%) contamination, we reduced the extraction radius 
to 6.5 pixels (15"). Such a region includes about 71% of 
the source PSF and the possible contamination from the 
nearby source within this region amounts to about 12% 
of its PSF. Given the faintness of the afterglow in this 
final part of the light-curve, this choice also improves the 
S/N ratio. 

The background level for the PC data was extracted 
in an annular region with an inner radius of 40 pixels 
and an outer radius of 150 pixels centered at the source 
position. To eliminate contributions from faint sources 
in the background region, we produced a 380 ks image 
by summing all of the observations, and searched this 
image for faint sources within the background annulus. 
In addition to the serendipitous source shown in Fig. ^ 
20 other sources were found with S/N ratio higher than 
3, all of which were located more than 50 pixels from the 
afterglow. The contributions from these sources were 
excluded from the background region. 
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Fig. 2.— BAT and XRT light curves of XRF 050416A. The XRT count rate (0.2 - 10 keV) was converted into flux units by applying a 
conversion factor derived from the spectral analysis. The solid line represents the best-fit model with the doubly-broken power law. The 
dashed line is the extrapolation of the XRT best fit model prior to the first XRT observation. The BAT light curve was extrapolated into 
the 0.2 — 10 keV XRT energy ba nd by converting th e 15 — 150 keV BAT count rate through the single power law best fit model with photon 
index F = 3.1 presented in iSakamoto et all 12005d) . 



Data were binned in order to have a S/N ratio higher 
than 3. The source was not detectable after observa- 
tion # 13; data from observations 14-19 were summed 
together and provide a single 3 a upper limit value of 
2.3 X 10~^ count s~^. Figure |21 shows the background- 
subtracted light curve in the 0.2 — 10 keV energy band. 
The source is clearly fading with time. 

The XRT light curve decay is not consistent with a sin- 
gle power law (x?ed = 1-43, with 80 de grees of freedom, 
d.o.f.). A broken power law: 

^ > \ K T^B-QA^-OB for t > Tb 



(where and ae are the power-law slopes before and 
after the break time Tbj respectively) improves the fit, 
giving Xred ~ ^-"^^ C^^ d.o.f.). However, the residuals 
show a systematic trend. Adding a second break to the 



model: 

r K for t < Tb,i 

F{t) = \ K T^^f-"* t-"^ for Tb,i <t < n,2 

[ K Tb"f""* r^"c-aB ^-ac for ; > 2 

the fit further improves yielding Xrcd = 0-^1 ('''^ d.o.f.). 
The F-test for this model vs. the simple broken power- 
law gives a chance probability of 1.2 x 10~^. This last 
model reveals the presence of two breaks at 172 ± 36 s 
and at (1.450 ± 0.013) xlO^ s. Table shows the best 
fit results obtained with the three models. We also tried 
to fit the light curve with a single or broken power law 
allowing the reference time to be a free parameter. 
For the simple power law and the single break broken 
power-law Xred improved to 1.40 (79 d.o.f.) and 1.02 (77 
d.o.f.), respectively, but the best-fit model is not able to 
account for the initial part of the light curve decay and 
the resulting tg is before the burst onset [t^ — — 201 ± 
156 s and to = — 153±34 s, respectively). For the doubly- 
broken power law the improvement is not significant and 
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TABLE 2 

XRF 050416A light curve best fit parameters 



Parameter Single PL Broken PL Doubly-broken PL 



«A 


0.82 ± 0.02 


3.28 ± 0.04 


2.4 ±0.5 


Tb,i (s) 




103 ± 9 


172 ± 36 


Ob 




0.81 ± 0.02 


0.44 ± 0.13 


Tb,2 (10=^ s) 






1.450 ± 0.013 


«c 






0.88 ± 0.02 


xLi (d.o.f.) 


1.43 (80) 


1.20 (78) 


0.81 (76) 



Note. — ua, aB, and ac are the decay s lope s 
for the distinct phases of the hght curve (see S I3.2|l . 
Tb.i and rb,2 are the epochs at which the decay slope 
changes, measured from the XRF onset. The IM and 
LR points have been included in the fits. 



the value of is marginally consistent with the burst 
trigger time. 

In Fig. [2 we plot the data as well as the best fit model 
obtained with the doubly-broken power law, with = 0. 
We also show the model extrapolation back to the time of 
the trigger. Hereafter, we will refer to the time intervals 
t < Tb,i, Tb,i < t < rb,2, and t > Tb,2 as phases 'A', 'B', 
and 'C, respectively. 

Note that the second observation starts after the sec- 
ond break in the light curve (see fifth column in Table^), 
and all other observations contribute only to phase C. 
Observations from 3 to 13 correspond to the last seven 
points of the light curve shown in Fig.|51 and observations 
from 14 to 19 give the final upper limit. 

3.3. SPECTRAL ANALYSIS 

Detailed spectral analysis could be performed only for 
the LR data in settling mode acquired before the start of 
the first observation and on the PC data of the first and 
second observations; the following 11 observations did 
not add significant statistics to the spectra, thus they 
were not included in the analysis. The same holds for 
the WT data of the first observation consisting of only 
31 (background-subtracted) photons. 

Given the presence of three different phases in the 
XRF 050416A light curve decay, we first checked for spec- 
tral variability during the afterglow evolution. The LR 
spectrum and background were extracted as described 
in § 13.21 The two PC spectra for the phases A and B in 
Fig. 121 were extracted from the first observation with the 
same annular region used for the timing analysis, while 
the PC spectrum for the phase C was extracted from the 
remaining part of the first observation and the second 
observation (i.e. up to 2.8 x 10^ s) using the circular 
region of 30 pixel radius. The PC background spectra 
were extracted from the same region as for the timing 
analysis. 

We fitted the spectra from phases A, B, and C sepa- 
rately with absorbed power laws. For phase A, the LR 
and PC spectra were fitted together leaving the normal- 
ization parameters free to take into account the different 
rate level due to the light curve decay. The four spectra 
are shown in Fig. |21 together with their absorbed power 
law best fits. 

The best fit results show evidence for spectral varia- 




0.2 0.5 1 2 5 10 

energy (keV) 



Fig. 3. — Phase A, B, and C spectra converted in the E^N{E) 
(or uFij) representation together with their best fit absorbed power- 
law models. The panels showing residuals refer to phase A, B, and 
C, respectively from top to bottom. 

tion among phases: the emission in phase A is signif- 
icantly softer than in the phases B and C. In the lat- 
ter two phases the best fit photon indices were consis- 
tent within the errors. Spectra extracted from phases B 
and C were, therefore, summed together and their an- 
cillary response files, obtained from different extraction 
regions, were weighted according to the relative exposure. 
The fit of the resulting spectrum gave an absorption col- 
umn density of {2.6~!_o'^) ^ ^^^^ cm~^ and a photon in- 
dex r = 2Mtom, with a reduced of 1.2 (83 d.o.f.). 
The observed column density is significantly larger than 
the Galactic value (0.21 x 10^^ cm~^). We therefore 
checked for intrinsic absorption in the host galaxy by 
adding a redshifted absorption component (zwabs model 
in XSPEC yl 1-3.1) with the redshift fixed to 0.6535 
ijCenko et al.l I2005q) and the Galactic absorption col- 
umn fixed to 0.21 x 10^^ cm~^. The fit gave a value 
cm for the additional column den- 
sity with an improvement of the reduced fo 1.0 (83 
d.o.f.). The phase A spectra, that also showed a column 
density significantly higher than the Galactic value, were 
fitted again with the addition of a redshifted absorption 
component. Because of the low statistics of the phase 
A data, the column density was kept fixed to the value 
obtained from the phase B-l-C spectral fit and the best 
fit photon index was F = S.OtoA (xL = 0-6; 5 d.o.f.). 
Table O shows the final results of the spectral analysis. 

4. DISCUSSION 
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TABLE 3 

XRF 050416A XRT spectral fit results 



Parameter Phase A Phase B+C 



Galactic column density (lO^^ cm ^) 


0.21 (frozen) 


0.21 (frozen) 


Host column density (lO'^^ cm~-^) 


6.8 (frozen) 


6.8tl:0 


Photon index 


O Q + 0.3 

•^■"-0.4 




N (photons keV~i cm"^ s~i at 1 keV) 


0.14 ±0.03 


(2.7 ±0.3) X 10"* 


0.2 - 10 keV flux (erg cm-^ s"!) 


(1.9 ±0.1) X 10-10 


(1.7 ±0.2) X 10-12 


0.2 — 10 keV luminosity (erg s"^) 


5.7 X 10*^ 


3.1 X lO*'^ 


X?ed (d-O.f.) 


0.6 (5) 


1.0 (83) 



Note. — The intrinsic column density for the LR data was 
held fixed to the best fit value found from the PC spectrum. The 
(isotropic) luminosity was calculated for a redshift z = 0.6535, with 
Ho=70 km s~^ Mpc~^, nm=0.3, and Qa=0-7- Unabsorbed fluxes and 
luminosities reported for the PC data are averaged over long time in- 
tervals: accurate instantaneous values for the unabsorbed flux can be 
derived from Fig. |21 



We have presented a detailed analysis of the X-ray 
afterglow of XRF 050416A. The prompt emission of 
this burst lasts ~ 2.5 s and is characterized by a first 
peak followed by a second much weaker one. The 
event be longs to the short tail o f the long GRB pop- 
ulation ijKouveliotou et all Il993j) . It could also be 
consistent with t he th ird class of bursts identified by 
iMukheriee et alJ l)1998D through multivariate analysis on 
the BATSE catalog. This class consists of bursts of 
intermediate duration and fluence as compared to the 
standard classes of short/faint/hard and long/bright/soft 
GRBs. The average energy distribution of the prompt 
emission is soft (well fitted by a power law with /3 = 
r - 1 = 2.1 ± 0.2 at the 90% confidence level or by a 
Band model with OfBand = —1, Ep = 15.612 7 ^^"^j and 
/^Band < —3.4 at the 68% confidence level) with signifi- 
cant evidence of hard to soft evolution within each peak 
(jSakamoto et aljr2005cj) . The gamma-ray spectral dis- 
tribution of this burst, with fluence in the X-ray energy 
band 2 — 30 keV larger than the fluen ce in 30 — 400 keV 
band, classifies it as a n X-ray flash ijLamb et alJl2005t 
ISa,ka,mntn et aTll2005l). 

XRT monitored the XRF 050416A X-ray emission from 
~ 64.5 s up to 74 d after the BAT trigger. XRF after- 
glows have been rarely de tected in the past (XRF 011030, 



XRF 
XRF 



020427: 
030723: 



Bloom eta] 200.-^ 
Butler et al.| I200F 



Levan et al. 200^ 
XRF 040701: Wo^ 



2(M XRF 050215B: iLevan et al.[l2?)05k XRF 050315 



Vauehan etaJ] I2006t XRF 050406: IBurrows et al 



2005bHRomano et a l."2006a': XRF 050824~'K rrmm et al 
2no5r The exceptionally long observational campaign 



of XRF 050416A has provided us with a unique data set 
and allowed one of the most accurate spectral and timing 
analyses ever performed for an XRF afterglow. 

The XRF 050416A hght curve of the first 10^ s af- 
ter the trigger (Fig. ^ is fairly smooth and similar in 
shape to other Swi ft-detected XRF and GRB afterglows 
like XRF 050315 liVauehan et alJl20(p ) or GRB 050319 
lClia!2maiiQ_£LaD |2QQ 6al; see also iN ousek et al. 2006; 
iChincarini et alJl2005tio'Brien et alJl2006.1 . It shows ev- 
iden ce of three different phases (A, B, and C according 
to § 13. 2() , each of them characterized by a distinct de- 



cay slope (see Fig. |21and Table EJ- At the beginning of 
the XRT observation the light curve shows a steep decay 
(aA ~ 2.4), followed by a short flat phase (ae ~ 0.44) 
and then by a third long- lasting phase with a more rapid 
intensity decline (ac ~ 0.9). We also found that the late 
extrapolation of the phase C decay is consistent with the 
flux upper limit measured 65 — 74 d after the prompt 
emission (observations # 14 — 19). There is no evidence 
of X-r ay flares as seen i n XRF 050406 llRomano et al 
12006a), GRB 050502B ( Burrows et al. '2005bt iFalcom 
12006), GRB 050607 ( Pagani et al. 200 6), and many other 
events. 

The XRT spectra show significant excess absorption in 
the rest frame of XRF 050416A (A^h ~ 6.8 x lO^i cm-^) 
and an energy distribution significantly softer in phase 
A {/3 — 2.0 ± 0.4) than in phases B and C, which have 
consistent spectral slope with a weighted average energy 
index /3 = 1.04 ± 0.05. This may indicate a different 
emission process acting during the initial phase of the 
XRF 050416A hght curve. 

4.1. PHASE A 

In the internal/external shock scenario 
ijRees fc Meszarosi .1994) . the tails of GRB peaks 
are expected to be cau sed by the "hi gh lat itude effect" 
l|Kumar fc PanaitesciJ Hmi iDermeJ I200I . If a rela- 
tivistic shell of matter suddenly stops shining, a distant 
observer receives photons emitted from increasing off- 
axis angles at later times due to the longer travel path. 
When the observed frequency is above the synchrotron 
cooling frequency (which is usually the case in the 
X-ray band), then the decay index expected for the 
observed light curve is a = 2 -I- /3, where /3 is the spectral 
index measured during the decay. The decay slope of 
phase A (aA = 2.4 ± 0.5) is definitely lower than the 
value expected for high-latitude emission (^ 4.0 ± 0.4). 
However, the high-latitude effect only provides an upper 
limit to the decay slope, since it assumes that the shell 
emission stops abruptly after the initial pulse. Residual 
emission may still be present from the shocked shells, so 
that slower decline rates are possible. The high-latitude 
emission would in this case contribute only a small 
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fraction of the overall radiation. Another possibility is 
that the X-ray band was below the cooling frequency: 
in this case the expected decay slope of high latitude 
radiation would naturally be shallower than 2+f3. 

In this scenario the first break in the X-ray light curve 
would be due to the emergence of the afterglow light af- 
ter fading of the prompt emission (see Sect. 14. 2|) . The 
forward shock may therefore contaminate the tail emis- 
sion. This component would however make the spectrum 
harder. Since most of the phase A counts come from the 
LR data, that correspond to the very first point in the X- 
ray light curve, we expect the afterglow to provide little 
contribution 10%) at this time. 

4.2. PHASES B AND C 

According to iZhang et al.l l)2006(l . the standard inter- 
pretation of the flat decay slope during phase B and 
of the second tem poral break in the afterglow involves 
refreshed shocks l|Sari fc Mcszaros 2000). In the ini- 
tial stages of the fireball evolution the forward shock, 
which emission produces the X-ray afterglow, may be 
continuously refreshed with the injection of additional 
energy. This could happen either because the cen- 
tral engine still emits continuous ly ijDai fc Lul 119981: 
IZhang fc Meszarosll200lt iDail I2004D . or because slower 
shells emitted at the burst time c atch up the fire- 
ball w hich has already decele r ated llRees fc MeszarosI 
19981 IPm aitescu et all 119981 iKumar fc Pira nI '200rt 
Sari fc Mes zaros 200n HZhang fc M eszaros 20021. Within 
this scenario, a flat decay of the afterglow is expected as 
the refreshed forward shock decelerates less rapidly than 
in the standard case. A transition to the standard after- 
glow evolution (i.e. a break) with no remarkable spec- 
tral changes is also expected when the add i tional energy 
supply ends. According to the'Zhan g et al.l 1)200(1) analy- 
sis, the various refreshing mechanisms can be character- 
ized by an effective index q < 1 such that a decay slope 
ttinj = (1 -|- q/2)(j -|- g — 1 is expected for the afterglow 
Hght curve until injection stops. For XRF 050416A, with 
(3 ~ 1.0 and ainj = ckb ~ 0.44, we can derive q ~ 0.3. 

Note that the phase C decay slope and spectral in- 
dex are marginally consistent with ac ~ (3p — 2)/4 and 
(3c = pI'^ for p ~ 2. This is expected for a fireball prop- 
agating in a uniform interstellar medium or in a stel- 
lar wind environment when the synchrotron cooling fre- 
quency is b elow the X-ray region, and before the jet break 
jSari et all [1993 (Chevaher fc Lil 120001) . Since phase C 
remarkably continues uninterrupted until the last XRT 
detection 42 days after the burst, this interpretation im- 
plies the absence of both a cooling break (expected in 
case of wind environment) or a jet break in the X-ray 
afterglow. Wind and magnetic field parameters can be 
easily adjusted to delay the cooling break after the end 
of the observational campaign of XRF 050416A, thus the 
circumburst environment cannot be distinguished based 
on X-ray data alone. 

4.3. JET MODELS 

The long and well sampled X-ray light curve of 
XRF 050416A is one of the best suited for testing more 
realistic jet models. The leading jet models for XRFs in- 
clude an on-beam unifor m jet with a very wide opening 
angle l|Lamb et alJl2005j) . an off-beam uniform jet with 
the line of sight outside the jet edge ijYamazaki et al.l 



12003(1 ■ a structured Gaussian-like jet with the viewing an - 
gle outside the bright Gaussian core l)Zhang et alJ2004'ah . 
and a two-component jet with the line of sight on the 
less energetic wider beam such as the co coon surround- 
ing a coUapsar jet ((Zhang et al.ll2004bD . The segment 
C displays a "normal" afterglow decay without signifi- 
cant features. This rules out the off-beam uniform jet 
which predicts an i nitial fast rise and a rapid decay (e.g. 
iGranot et al.ll2?)0 3). The lack of a rebrightening feature 
also greatly constrains the two-component jet. We there- 
fore only focus on two possibilities, i.e. an on-beam uni- 
form jet and an off-beam structured jet. 

In this section we describe preliminary results of our 
efforts to model the light curve in the phase C through 
numerical simulations. We model the synchrotron radia- 
tion powered by the external shock propagating in a uni- 
form medium with proton number density n. Our three 
dimensional code, which was developed by J. Dyks, can 
calculate afterglow emission of an arbitrary axially sym- 
metric jet observed at any viewing angle. The code has 
been used to mod el the curvature effect in structured jets 
l(Dvks et al.l2005() and more generally the GRB afterglow 
light curves (described in Zhang et al. 2005). The code 
takes into account all kinematic effects that affect the ob- 
served flux (e.g. Doppler boost, propagation time delays) 
in the way described by Salmonson (2003). The radial 
dynamics of the outflow is followed well i nto the non- 
relativistic regime using the equations of iHuang et alJ 
(2000). Sideways expansion is neglected. The evolu- 
tion of the electron energy spectrum is followed with the 
simplified analytical method (Fan et al. 2004), except for 
the top-hat (i.e. uniform jet) case when the low run time 
has allowed exact integration of the continuity equation 
l(Moderski. Sikora fc Bulikl2000(l . In the calculations dis- 
cussed below, we have used the exact spectrally derived 
value of p = 2/3c = 2.08, which assumes the cooling fre- 
quency was below the X-ray band. 

Light curves calculated for a uniform jet viewed on- 
beam are notably steeper than the data (the predicted 
slope is a = 1.1, to be compared with the observed value 
UQ = 0.88 ± 0.02). The one shown in Fig. H (solid fine) 
has been calculated for the half opening angle of the uni- 
form jet = 17°, the viewing angle C = 0, the total 
explosion energy of the two-sided outflow E ~ 2 x 10^^ 
erg, the bulk Lorentz factor Fbik = 80, the proton num- 
ber density of the external medium n = 0.1 cm^"^, and 
the electron and magnetic field energy equipartition pa- 
rameters, ee and eb, equal to 0.1 and 0.01, respectively. 
A "contrived" version of structured jets, with uniform 
distribution of Fbik as a function of 9, but with the ex- 
plosion energy e{0) increasing towards the jet axis can 
produce light curves that are less inconsistent with the 
data (dashed line in Fig. ^J. This is thanks to the well 
known flattening that appears near the jet break time 
for large viewing angles C > 6*0010, where ^coro is the 
half opening angle of the most energetic central part of 
the outflow. However, the light curves typically have a 
slightly concave shape that is not observed in phase C. 
The dashed line in Fig. ^ is for a structured jet with the 
uniform Fbik, and a Gaussian profile of e{6) with stan- 
dard deviation a = ^corc = 4° and ( — 10°. The other 
parameters are the same as in the top-hat case. 

Standard theoretical parameter values of Fbik =100, 



The Afterglow of XRF 050416A 



10-9^ i 



1 1 1 1 1 — I — I I 1 1 1 1 1 



I 1 1 1 1 1 1 — I — I I 1 1 1 1 1 1 — I — I I 1 1 1 1 1 




10^ 10^ 
Time (seconds from trigger) 



Fig. 4. — Two model light curves plotted on the XRT data. The solid line shows the results for a uniform jet viewed on-axis = 0°) 
with opening angle djgt = 17°, FbU; = 80, and E = 2 X 10^^ erg. We have assumed n = 0.1 cm~^, and p = 2.08. The dashed line shows 
the results for a structured jet with a Gaussian distribution of outflow energy and with a uniform distribution of Fbiij, viewed off-axis 
((^ = 10°). The la width of the Gaussian was 4°. The other parameters are the same as in the uniform case. 



n — 1 cm^^ and an assumption oi E = 2 x 10'^^ erg 
result in models in which the cooling frequency is be- 
low the X-ray band throughout phase B, with the result 
that the model slope is steeper than the observed slope 
in that phase. The flatter slope of phase B could be in- 
terpreted as a "cooling break" if the cooling frequency 
begins above the XRT energy band and crosses it at the 
phase B / phase C transition. This would require lower 
Tbik ~ 70 — 80, lower n ~ 0.1 cm~'^, or larger explo- 
sion energies E > 10^^ erg. The low value of n and 
large E could also help explain the late jet break time 
without requiring a large solid angle for the uniform out- 
flow. However, the presence of a cooling break at the 
phase B / phase C transition should be associated with 
a strong spectral change, which is not observed. In fact, 
the modeled energy spectral index in this case changes 
by 0.5 at this transition (see Fig. O, in contrast to the 
observations. The hard model spectrum in phase B could 
in principle be compensated for a contribution from the 
softer, prompt emission, but this seems contrived, as it 
would have to balance perfectly to produce no change 



in the observed spectral index. The interpretation given 
above for phase B in terms of refreshed shocks seems 
more plausible. 

5. GONCLUSIONS 

The steep initial phase of the X-ray light curve of 
XRF 050416A can be interpreted as the tail of the 
prompt emission. The rest of the light curve shows evi- 
dence of forward shock refreshing mechanisms acting up 
to about 1450 s since the trigger and then an uninter- 
rupted decay with no signature of jet breaks up to at 
least 42 days. 

The overall phenomenology of XRF 050416A and other 
Swift detected XRFs provides evidence that both GRBs 
and XRFs arise from th e same pheno menon l ^^rnb ^1^ all 
l2f)05l ISa,ka,moto etHI L2005d: D'Al essio PirnI 120051) . 
but the characteristics of XRF 050416A deviate some- 
what from model expectations. Our preliminary nu- 
merical simulations confirm that the late-time decay of 
XRT 050416A is slightly flatter than expected for on-axis 
viewing of a uniform jet. Modeling of off-axis viewing of 
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Fig. 5. — The temporal evolution of the energy spectral index /3 
calculated for models with n = 0.1 cm~^, E = 2 X 10^^ erg and 
three different values of the bulk Lorentz factor Fj^uj = 100, 80, 70 
(dashed, dotted, and solid). The index was calculated as an arith- 
metic average of 5 values uniformly distributed within the XRT 
energy band. 



Swift XRT results are showing that jet breaks ar e rare 
for b oth XRFs and G RBs (XRF 050315 . FVaughan et all 
"2006^ GRB 50318, ' Perri et alJ GRB 050505, 

HurketTeTaOBoOe; GRB 050525AlBlustin et al. 2005). 
Because of its relative slow decay rate and its excep- 
tionally long X-ray light curve, XRF 050416A provides 
a particularly stringent case for understanding jet col- 
limation and structure in GRBs. A jet opening angle 
^jetj^ 28° can be estimated through Eq. 1 in Sari et aO 
l)1999|) assuming ijot ^ 42 d (i.e. the time of the last 
detection in the XRT light curve), an isotropic energy 
Ei^o r$ 1.2 X 10^2 gj.g (obtained from the BAT power law 
fit) and standard values of the proton number density of 
the external medium n — i and the radiative conversion 
efficiency of the burst rj^ = 0.2. This opening angle far 
exceeds the jet angle of 5° — 10° generally considered to 
be typical for GRBs. Whether this is related to the soft 
nature of this event requires a larger sample of jet breaks 
for both GRBs and XRFs. 



a structured jet with a uniform bulk Lorentz factor and 
a Gaussian internal energy distribution also could not 
reproduce the X-ray light curve satisfactorily for reason- 
able values of the model parameters. More realistic af- 
terglow models should be investigated in detail to under- 
stand these effects and to fully explain the broad range of 
peak energies observed from XRFs through X-Ray Rich 
bursts to classical GRBs. 
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APPENDIX 
XRF 050416A AND THE BAND MODEL 

O'Brien et al. (2005) have shown that for most GRBs the X-ray afterglow light curves smoothly connect with the 
end of the prompt emission. To show this, it is necessary to compute the expected flux in the XRT energy range 
due to the prompt emission. In Fig. [3 we show that the the backward extrapolation of the X-ray light curve of 
XRF 050416A to the time of the trigger smoothly joins with the extrapolation of the BAT light curve in the XRT 
energy band obtained extending the BAT best fit power law with F = 3.1 down to soft X-ray energies. However, the 
XRF nature of the source implies that the BAT light curve in Fig. [21 represents only an upper limit to X-ray emission 
du ring the prompt phase. A BAT light curve extrapolation computed by adopting the Band best fit model presented 
bv lSakamoto et alJpOOS^) . or any Band model with —2 < ctBand < —0.5 and Ep 16 keV, predicts a flux significantly 
below the backwards extrapolation of the phase A XRT light curve (more than one order of magnitude) . 

This would suggest that the steep decline phase is not tied to the prompt emission, but is, perhaps, due to a further 
burst peak or a flare (too faint to be detected by BAT) which occurred during the data gap not covered by XRT. 

It is however interesting to note that the spectral index of the phase A XRT light curve is very soft, and is consistent, 
within the errors, with that of the high- energy slope of the main burst spectrum. The phase A light curve could therefore 
be indeed the prompt emission tail (as with most Swift GRBs), provided that the peak energy had shifted red- ward 
of the XRT range by the beginning of the XRT observation. Indeed, softening of GRB spectra is a common property 
of their prompt emission (Ford et al. 1995; Romano et al. 2006b). Ghirlanda et al. (2002) showed that in some cases 
this behavior is due to the lowering of Ep. Thanks to Swift, this phenomenon might have now been observed over a 
much wider temporal and spectral span. 

The XRF 050416A BAT spectra sho wed spectral evolution l|Sakamoto et al.ll2005cD . with the hardness ratio decreas- 
ing during each burst peak. Fig. IA6I shows the expected values of the (25 — 50 keV)/(15 — 25 keV) hardness ratio 
in simulated BAT observations of a source with a Band spectral energy distribution located at the same detector 
position as XRF 050416A. The hardness ratio is plott ed as a function of E-p for different choices of the indices ctBand 
and /^Band, covering their typical ranges of variation l)Preece et al J 1200(1(1 . Simulations of BAT spectra with a Band 
shape with aeand = —1, /3Band — —3, and Ep varying from 0.1 to 50 keV show that the observed hardness ratio of 
the XRF 050416A BAT peaks (which decreased from ~ 1 to 0.4) c an be reproduc ed by Ep evolving from ~ 30 
to ^ 10 keV, consistent with the average value Ep ^ 16 keV (Sakamo to et al.ll2005d) . A similar behavior could be 
obtained also allowing an evolution of aeand (from aeand > — 1 to aeand < — l)j while pure /3Band evolution (keeping 
fixed Ep — 16 keV and ctBand = —1) cannot easily account for the observed hardness ratio range. 

To summarize, the BAT and XRT data are consistent with the peak energy evolving from 30 keV (during the 
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Fig. A6. — Plots of (25 — 50 keV)/(15 — 25 keV) hardness ratio values from simulated BAT observations of a source with Band spectral 
energy distribution at the same position on the detector array as XRF 050416A. Top: Hardness ratio as a function of Ep, for different 
values of /^Band with OBand fixed to —1. Bottom: Hardness ratio as a function of Ep, for different values of aeand with /Ssand fixed to —3. 



second GRB peak) down to 1 keV (at the beginning of the XRT observation) . This would roughly correspond to 
Ep oc t~^. Indeed, fixing the intrinsic absorption to the value obtained by the fit of the Phase B and phase C spectra, 
the phase A XRT spectrum may be fit by a Band model with Ep S 1 keV (although the data do not require this). 
Otherwise, Ep could be even lower, close to or below the XRT range (so that the Ep evolution would be faster). 
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